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Abstract 

Increasing integration density of electronic components 
has exacerbated the thermal management challenges facing 
electronic system developers. The high power, heat flux, and 
volumetric heat generation of emerging devices are driving 
the transition from remote cooling, which relies on conduction 
and spreading, to embedded cooling, which facilitates direct 
contact between the heat-generating device and coolant flow. 
Microgap coolers employ the forced flow of dielectric fluids 
undergoing phase change in a heated channel between 
devices. While two-phase microcoolers are used routinely in 
ground-based systems, the lack of acceptable models and 
correlations for microgravity operation has limited their use 
for spacecraft thermal management. Previous research has 
revealed that gravitational acceleration plays a diminishing 
role as the channel diameter shrinks, but there is considerable 
variation among the proposed gravity-insensitive channel 
dimensions and minimal research on rectangular ducts. 
Reliable criteria for achieving gravity-insensitive flow boiling 
performance would enable spaceflight systems to exploit this 
powerful thermal management technique and reduce devel- 
opment time and costs through reliance on ground-based 
testing. In the present effort, the authors have studied the 
effect of evaporator orientation on flow boiling performance 
of HFE7100 in a 218 um tall by 13.0 mm wide microgap 
cooler. Similar heat transfer coefficients and critical heat flux 
were achieved across five evaporator orientations, indicating 
that the effect of gravity was negligible. 


Introduction and Motivation 

Increasing functionality and miniaturization of modern 
and emerging electronic components has exposed the 
limitations of the current remote cooling paradigm, which 
relies on conduction and spreading across multiple interfaces 
to dissipate waste heat. The large temperature gradient 
between the heat source and sink that results from remote 
cooling has led to electronic systems that are thermally 
limited, operating below the electrical capability of the device 
technology they exploit [1]. Embedded cooling overcomes 
these limitations by facilitating direct contact between the 
heat-generating device and coolant flow. Systems that enable 
the forced coolant flow to undergo phase change within the 
embedded channels provide additional benefits, such as higher 
heat transfer coefficients, smaller temperature gradients, hot 
spot mitigation, and lower pumping power requirements. 

Space missions include an array of electronic and power 
systems, many of which could benefit from embedded 
cooling. These systems include power electronics, lidar and 
radar systems, power generation systems (e.g., Rankine cycle 
power plants, radioisotope thermoelectric generators, and fuel 
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cells), and three-dimensional integrated circuits. The latest 
NASA Technology Roadmap [2] lists needs for (1) acquisi- 
tion and removal of heat fluxes greater than 100 W/cm’ over 
relatively small areas with tight temperature control (+ 1°C); 
(2) enhanced heat transfer surfaces with micro- and nano- 
scale features to enhance two-phase heat transfer at higher 
heat fluxes; and (3) high-capacity, two-phase heat transport 
systems for thermal control of large heat loads, such as those 
required by Rankine cycle power plants. 

A key benefit of embedded two-phase cooling for space 
missions is the ability to deliver waste heat from the heat 
source to the radiator with little temperature drop. This 
enables the radiator to operate at higher temperatures, 
reducing the required heat rejection area and thereby reducing 
the radiator size and mass. For the same heat load in low earth 
orbit, a radiator operating at 20°C must be 34% larger than 
one operating at 40°C. Other benefits of two-phase mechani- 
cally pumped fluid loops include longer transport distances 
with potential use of multiple evaporators and condensers and 
precise flow rate control, which increases heat flux limits and 
enables shutdown of the cooling system when necessary (e.g., 
when the heat generating device is unpowered or the space- 
craft enters survival mode) [3]. 

One of the significant barriers to the widespread use of 
two-phase microcoolers is the complex nature of flow boiling, 
particularly for microgravity applications for which only 
limited experimental data is available. The complexity can be 
reduced through the use of a single, low aspect ratio channel 
(1.e., a microgap) rather than an array of parallel 
microchannels. This configuration mitigates flow instabilities 
and possible flow reversals as the vapor can expand both 
spanwise and downstream [4]. However, two-phase microgap 
coolers are not widely employed for spacecraft thermal 
management due to the lack of data and correlations for 
microgravity environments and the absence of acceptable 
models that could enable extrapolation of heat transfer and 
flow behavior from available terrestrial data. 

Application of physics-based flow regime maps and an 
insightful reading of the two-phase literature reveals that 
gravitational acceleration plays a diminishing role as the 
channel diameter shrinks, thus facilitating more reliable 
extrapolation from existing databases. The goal of the present 
effort is to characterize the fluid physics governing two-phase 
flows in heated miniature and microscale rectangular ducts, 
with emphasis on methods for minimizing the effect of 
gravity in such flows. Experimental validation of gravity- 
insensitive behavior would enable spaceflight systems to 
exploit this powerful thermal management technique and 
reduce development time and costs through reliance on 
ground-based testing. 


Microscale Definitions 

Microscale in the two-phase flow sense refers to flows for 
which the controlling mechanisms differ from those at the 
macroscale, with the influence of surface tension and shear 
forces increasing and that of gravity diminishing [5]. Accord- 
ingly, a literature review has been performed to find criteria 
for the transition to microscale two-phase flow behavior, 
beyond which gravity effects should be negligible. Many 
microchannel definitions have been proposed, including those 
based on channel geometry alone and those relating to a 
combination of channel size and fluid properties such as 
surface tension and the density difference between the liquid 
and vapor. Recent definitions have incorporated additional 
terms, such as viscosity, mass flux, and contact angle. 


Definitions based on Geometry Alone 

Mehendale, Jacobi, and Shah [6] classified hydraulic di- 
ameters greater than 6 mm as macroscale, those from | to 6 
mm as compact, those from 0.1 to 1 mm as mesoscale, and 
those from 0.001 to 0.1 mm as microscale. Kandlikar [7], 
based on manufacturing techniques, defined conventional 
channels as those with hydraulic diameters less than 3 mm, 
minichannels as those with hydraulic diameters less than 3 
mm and greater than 200 um, and microchannels as those 
with hydraulic diameters less than 200 um. Such definitions 
based on geometry alone capture the compact nature of the 
channels, but fail to characterize the underlying change in the 
two-phase fluid behavior. 


Definitions based on the Edtvos Number 

Many microchannel definitions are some variation of the 
ratio between the gravitational forces (Ap - g - D,,) and surface 
tension forces (o/D;,) within the system. A useful way to 
assemble and compare these criteria is to reduce each one to 
its equivalent Edtvés number, Eo, defined as 


. . 2 
= Ap g Dy (1) 
oO 


Eo 


where Ap is the density difference between the liquid and 
vapor phases of the working fluid, g is the acceleration due to 
gravity, D;, is the hydraulic diameter (Dy, = 4: A,/P where 
Ag is the cross sectional area and P is the wetted perimeter), 
and o is the surface tension of the working fluid. Low values 
of the Edtvés number indicate that surface tension dominates. 

In 1963, Suo and Griffith [8] studied mixtures of liquid 
water, n-heptane, and n-octane with gaseous air, nitrogen, and 
helium in tubes with radii from 0.514 to 0.795 mm. They 
observed that liquid slugs in adiabatic horizontal tubes 
exhibited an insensitivity to gravity when 
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where {2 represents the ratio of gravitational forces to surface 
tension forces, p, and pg are the liquid and vapor densities, 
respectively, r is the tube radius, and 4, and yy are the liquid 
and vapor dynamic viscosities, respectively. By replacing the 
liquid density with the density difference between the liquid 
and vapor phases, which introduces only a small error due to 
the constraint of p,;/p, >> 1, and using the tube diameter 
rather than the radius, the Suo and Griffith criterion for the 
microscale transition can be rewritten as Eo < 0.88. 

In 1992, Brauner and Moalem Maron [9], using the strati- 
fied/nonstratified transitional boundary as their guide, 
identified microchannels as those with Eo < 42*. Brauner 
[10] had previously provided the same definition to define 
channels subjected to Earth's gravity that exhibited similar 
behavior to that expected in microgravity. 

In 1997, Kew and Cornwell [11] found that established 
correlations predicted heat transfer coefficients reasonably 
well for flow boiling of R141b in tubes with diameters of 3.69 
and 2.87 mm, but performed poorly when applied to a tube 
with a diameter of 1.39 mm. They suggested that microscale 
behavior was exhibited when Co = 0.5. Cornwell and Kew 
[12] had previously introduced the confinement number, Co. 
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The E6tvés and confinement numbers are therefore related 
through the following expression. 


Eo =— > (6) 


Thus, the Kew and Cornwell criterion for microchannels can 
be rewritten as Eo < 4. 

In 1999, Triplett et al. [13,14] studied flow patterns, void 
fraction, and pressure drop of air-water mixtures in circular 
channels with diameters of 1.10 and 1.45 mm and semi- 
triangular channels with hydraulic diameters of 1.09 and 1.49 
mm. They suggested microchannels are those with hydraulic 
diameters smaller than the Laplace constant, La. 
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The Triplett et al. criterion can be rewritten as Fo < 1. 

In 2003, Li and Wang [15] studied the effect of tube size 
on condensation in mini- and micro-tubes using an analytical 
model. They found that flow patterns in small channels were 
limited to annular, lengthened bubble, and bubbly flow. They 
argued that the effect of gravity on the flow regime could be 
ignored and the thickness of the annular liquid film around the 
tube would deviate by no more than 1% when 
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which can be rewritten as Eo < 0.5. Follow on work in 2006 
by Cheng and Wu [16] indicated that gravity-insensitive 
behavior was seen in systems with Eo < 0.05. The authors 
also found that surface tension becomes dominant and the 
gravitational effect is small in systems with 0.05 < Eo < 3. 

In 2007, Ullman and Brauner [17] classified channel sizes 
based on flow patterns and their transitions. They found that 
modified flow pattern transitions were required for systems 
with Eo < 1.6. The authors also noted that — in calculating the 
E6tvés number for noncircular tubes — the hydraulic diameter 
might not represent correctly the significance of gravity and 
surface tensions forces (a concern also voiced by Baldassari 
and Marengo [18]). For rectangular ducts configured with 
their width and axial length normal to the gravity vector, 
Ullman and Brauner proposed that 
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where H is the channel height and W is the channel width. 
This variation of the Eétvés number accounts for the gravity 
term scaling with the channel height (Ap-g-H) and the 
surface tension term scaling with the channel width (o/W). 
The effect of this modification for HFE7100 in a 13.0 mm 
wide channel for channel heights from 1 to 10,000 um is 
shown in Figure 1. HFE7100 fluid properties were obtained 
from Rausch et al. [19]. The horizontal black lines represent 
the maximum (477) and minimum (0.05) values of the 
microscale Eé6tvés numbers proposed by various authors. The 
area of interest (i.e., the shape confined by the red, blue, and 
black lines) covers a range of channel heights spanning four 
orders of magnitude, revealing the significance of using the 
appropriate formulation of the Eétvés number. To the authors' 
knowledge, no experimental studies have been performed to 
assess which formulation is appropriate for rectangular ducts. 
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Figure 1: Large variation in the Eétvés number results from 
differing formulations based on hydraulic diameter (Equation 1) 
and channel height and width (Equation 9) 


In 2011, Ong and Thome [20] studied flow boiling pat- 
tern transitions and compared the liquid film thickness 


between the top and bottom of the channel for tubes with 
diameters of 1.03, 2.20, and 3.04 mm and working fluids 
R134a, R236fa, and R245fa. They observed a gradual 
transition from macroscale to microscale. The authors 
proposed Co > 0.3 —0.4 as the boundary for macroscale 
behavior and Co > 1.0 as the upper boundary for symmetric 
microscale flow. These values correspond to Fo < 6.25 for 
the boundary for macroscale (using Co = 0.4) and Eo < 1 for 
the boundary for symmetric microscale flow. 

Figure 2 shows the threshold Eétvés numbers for 
microscale two-phase flow proposed by the various authors 
discussed in this section. The ratio between the maximum and 
minimum values is nearly 800. Some of the variation can be 
attributed to the variation among the parameters each author 
used to assess microscale behavior. Authors who studied the 
absence of stratified flow generally reported higher values and 
authors who studied the uniformity of the liquid film thick- 
ness in annular flow reported lower values. Nonetheless, there 
exists significant variation among the criteria even when 
authors report using the same phenomenon to determine their 
microscale criterion. The wide variation suggests that the 
Eétvés number alone may not capture adequately the 
transition to microscale two-phase flow. 
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Figure 2: Eétvés numbers for microscale behavior proposed by 
various authors span three orders of magnitude 


Definitions beyond the Eétvés Number 

Recent efforts have considered additional terms to assess 
the transition to microscale two-phase flow behavior. In 2010, 
Harirchian and Garimella [21] studied confinement effects in 
flow boiling of FC77 in 12.7 mm long channels with widths 
of 100 to 5850 um and depths of 100 to 400 um. The authors 
suggested that — in addition to channel size and fluid proper- 
ties — mass flux also governs bubble confinement. They 
observed confined flows when the convective confinement 
number, Bo°> - Re, was less than 160 with 
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where Bo is the Bond number, D = VA. and G is the mass 
flux. Earlier work by Harirchian and Garimella [22] revealed 
that the channel cross-sectional area played a critical role in 
determining microchannel heat transfer mechanisms and thus, 
the characteristic length used in their analysis is the square 
root of the channel cross section rather than the hydraulic 
diameter. The modification in calculating the Bond number 
produces the same value that Ullman and Brauner [17] 
proposed for rectangular ducts, assuming the channel is 
oriented horizontally (the Bond and Eétvés numbers have 
identical formulations in the context of the present study). The 
effect of the modification for HFE7100 at 101.3 kPa ina 13.0 
mm wide channel for channel heights from | to 10,000 um is 
shown in Figure 3. 
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Figure 3: Large variation in the convective confinement number 
results from formulations based on hydraulic diameter (Equa- 
tion 1) and channel cross section (Equation 10) 


In 2015, Tibirica and Ribatski [23], by combining a broad 
literature review with their own experiments using tubes with 
diameters from 1.00 to 2.32 mm, concluded that the key 
characteristics of microscale tubes are (1) the absence of 
stratified flows and (2) uniformity of the liquid film thickness 
along the tube perimeter for slug and annular flows in 
horizontal channels. Using an analytical model, they deter- 
mined the maximum diameter for which a liquid slug could 
exist in a tube under static conditions was 
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where @ is the contact angle and p, is the vapor density. Thus, 
stratified flow is not expected when Eo < 8- cos @. Based on 
their experiments and the model of Kandlikar [5], the authors 


indicated that uniform film thickness was obtained when the 
gravitational force is less than 5% of the surface tension force. 
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Thus, the uniform film thickness criterion proposed by 
Tibirica and Ribatski [23] matches the microchannel criterion 
of Cheng and Wu [16]. 


Summary of Microchannel Definitions 

The recent criteria for the transition to microscale behav- 
ior cannot be reduced to the Eétvés number alone. Therefore, 
to provide a comparison among all of the proposed microscale 
criteria, Figure 4 shows the range of diameters at which 
various authors have suggested the transition to microscale 
behavior begins for saturated water, R245fa, and HFE7100 at 
101.3 kPa (corresponding to saturation temperatures of 100.0, 
15.1, and 59.8 °C, respectively). The fluid properties for water 
and R245fa were obtained from REFPROP Version 9.1 [24] 
(which calculates properties for water according to Wagner 
and Pruss [25] and R245fa according to Lemmon and Span 
[26]) and those for HFE7100 were obtained from Rausch et 
al. [19]. The contact angles for water and R245fa are from 
[23] and for HFE7100 from [27]. It is assumed that the 
channels have circular cross sections, thereby providing 
identical values for the physical and hydraulic diameters. 

The most conservative estimates for the transition to 
microscale two-phase flow behavior in circular channels are 
0.336, 0.237, and 0.191 mm for water, R245fa, and HFE7100, 
respectively, and the most relaxed values are 15.7, 6.65, and 
5.36 mm, respectively. The large variation among the 
proposed criteria (spanning two orders of magnitude for 
circular ducts) suggests additional research is required to 
clarify the microscale transition and that there may be 
additional parameters that influence gravity-insensitivity. 

Some studies have been performed with the aim of better 
understanding the role of gravity on two-phase flow behavior 
(rather than parameters that provide gravity-insensitivity), but 
those efforts cannot be covered comprehensively in this paper 
due to length constraints. The authors refer readers searching 
for additional details to several review papers [18, 28-30]. 


Objectives 

The literature review revealed a wide range of criteria for 
predicting the transition to microchannel two-phase flow 
behavior. The large variation in the criteria, coupled with the 
unresolved issue of which length scales are appropriate for 
calculating the Edtvés for rectangular ducts, prevents the 
reliable use of existing criteria without further validation. 
Accordingly, a research program was initiated to: 

1. Characterize the parameters, including channel size, 
fluid velocity, fluid properties, flow regime, vapor quali- 
ty, and heat flux, that provide gravity-insensitive flow 
boiling performance in microgap coolers; 

2. Determine the appropriate metric for assessing gravity- 
insensitive flow boiling performance (e.g., flow regimes 
and their transitions; measurements such as heat transfer 
coefficients and pressure drop; or instabilities, such as 
critical heat flux and local dryout); and 


3. Assess the role of microgap channel length, in which the 
flow is never "fully developed," on two-phase flow 
characteristics and the onset of gravity-insensitive be- 
havior. 


As a first step in meeting these objectives, a test facility 
has been developed and preliminary experiments on the flow 
boiling performance of HFE7100 in a 218 um tall by 13.0 mm 
wide microgap have been performed. 
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Figure 4: Microscale diameters for circular channels for saturated water, R245fa, and HFE7100 at 101.3 kPa 
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An experimental apparatus was designed and assembled 
to provide two-phase flow test data with the evaporator in 
various orientations with respect to the gravity vector. The 
flow loop for orientation testing, a schematic of which is 
shown in Figure 5, supplies degassed fluid to the evaporator 
test section at the prescribed flow rate, temperature, and 
pressure. The working fluid is circulated by a gear pump with 
an electromagnetic drive (Micropump Series GA _ with 
EagleDrive). The preheater provides the desired liquid 
subcooling at the inlet to the evaporator test section. Fluid 
leaving the evaporator is condensed and/or subcooled via a 
plate heat exchanger (Lytron LL520G14) connected to a 
recirculating chiller (Julabo FP50). The circulator working 
fluid is distilled water. A 15 um sintered metal filter collects 
particulate contamination within the flow loop. The tubing 
material throughout the loop is stainless steel. A temperature- 
controlled, two-phase reservoir regulates the low-side 
pressure for the loop and compensates for the expansion and 
contraction of the working fluid during testing. 

The liquid flow rate is measured by a Pelton-type 
microturbine liquid flow sensor with a range of 20-200 
ml/min and accuracy of + 0.5% of the full scale rating 
(McMillan Flow Products Model 104). The fluid temperature 
is measured via four-wire, class 1/10 DIN, 100 Q resistance 
temperature detector probes at the preheater inlet, preheater 
outlet, condenser inlet, condenser outlet, and reservoir vapor 


space. The fluid temperature is measured via 0.51 mm 
diameter stainless steel type-T grounded thermocouple probes 
at the evaporator inlet and outlet. The absolute pressure is 
measured in the reservoir vapor space and evaporator inlet by 
silicon, strain gauge type transducers with ranges of 0 to 300 
kPa and accuracies of + 0.1% of the full scale rating (Honey- 
well FP2000). The differential pressure between the evapora- 
tor inlet and outlet manifolds is measured by a variable 
reluctance pressure sensor with replaceable pressure sensing 
diaphragms, each with an accuracy of + 0.25% of the full 
scale rating (Validyne Engineering DP 15). 

The working fluid is HFE7100, due to its saturation prop- 
erties (boiling point of 59.8 °C at 101.3 kPa), low freezing 
point (-135 °C), low electrical conductivity, non-toxicity, and 
non-flammability. The flow loop was designed to minimize 
the system leak rate due to concerns of air infiltrating the 
system, which could affect the fluid properties and boiling 
performance of the degassed fluid. Chen and Garimella [31] 
found that undegassed FC77, relative to degassed FC77, 
experienced higher pressure drops and more flow instability; 
"pseudo-boiling" (i.e., the formation of air-vapor bubbles that 
alter the flow field but do not enhance boiling heat transfer); 
and lower boiling incipience temperatures. Similar results 
were found by Muller-Steinhagen, Epstein, and Watkinson 
[32] and Sawada et al. [33]. Another compelling reason to 
ensure the degassed fluid remains free of air is that published 
properties are available exclusively for the degassed fluid. 
Weld, VCR, and ConFlat fittings and flanges were used 
wherever possible due to their very low leak rates. Prior to 


charging the loop, the working fluid was subjected to multiple 
freeze-pump-thaw cycles to remove non-condensable gases. 
The saturation pressure and temperature of the degassed fluid 
were measured and compared against reference data [34] to 
verify that the non-condensable gases had been removed. 
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Figure 5: Flow loop schematic 


The heat source for the evaporator is a 12.7 mm by 12.7 
mm by 0.6 mm silicon thermal test chip (TTC) mounted to a 
printed circuit board (Thermal Engineering Associates TTV- 
4102). The TTC provides uniform heating with a four-wire 
heater circuit design to eliminate parasitic heat losses in the 
supply wiring. The resistors cover more than 85% of the die 
area. Ten temperature-sensing diodes provide temperature 
measurements of the TTC. The diodes are located on the chip 
frontside. Measurements of the surface via confocal micros- 
copy revealed an average roughness of 0.031 um + 0.018 um. 

The diodes on the TTC are calibrated using a constant 
temperature air oven and a four-wire, class 1/10 DIN, 100 Q 
resistance temperature detector probe. The temperature of the 
oven was held constant until consecutive readings of tempera- 
ture and resistance, taken 10 minutes apart, varied by less than 
0.1 °C and 1 Q, respectively. 

The evaporator assembly, shown in Figure 6 and Figure 
7, is composed of a thermal isolator base (not shown), base 
plate, thermal test vehicle, fluid enclosure, polycarbonate 
cover, top plate, and thermal isolator cap (not shown). The 
fluid enclosure includes the fluid inlet and outlet taps and 
manifolds, pressure taps for absolute and differential pressure 
measurements, and the evaporator inlet and outlet fluid 
temperature probes. Replaceable, transparent polycarbonate 
covers are inset to the fluid enclosure. The microgap is 
located between the exposed face of the TTC and the interior 
face of the polycarbonate cover. Isolators minimize thermal 
losses to the ambient air and mounting plate. 

A Keyence VHX-5000 digital microscope with a 100 to 
1000x objective and precision stage was used to measure the 
microgap height. The height of the stage was adjusted in 10 


uum increments until the top of the silicon TTC was in focus. 
Using this height as the zero reference, the height of the stage 
was then adjusted until the bottom of the polycarbonate cover 
was in focus to determine the microgap height. The uncertain- 
ty of this measurement technique is approximately one full 
increment in the vertical direction due to the reliance on the 
operator to declare the image in focus (+ 10 pm). For the 
present assembly, nine measurements were taken, which 
produced an average microgap height of 229 um, standard 
deviation of 7 um, and range of 220 to 240 pm. 
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Figure 6: Evaporator assembly 


Figure 7: Axial cross section of evaporator assembly (to scale) 


The optical measurements of the microgap height were 
performed at 22 °C with the interior cavity of the flow 
enclosure open to atmospheric pressure. During two-phase 
testing, the flow enclosure and polycarbonate temperature was 
55 to 60 °C based on a nominal inlet saturation temperature of 
61 °C, liquid inlet subcooling of 2 to 6 °C, and convective 
losses to the ambient air. The pressure effect on the microgap 
height was neglected, as the saturation pressure of HFE7100 
at 60 °C is very close to atmospheric pressure. Accounting for 
the thermal expansion of the polycarbonate cover island and 
stainless steel enclosure up to the ledge on which the polycar- 
bonate cover rests and assuming a temperature rise of 35 °C 


results in an estimated reduction in the microgap height of 11 
um. Thus, the predicted microgap height is 218 pm. 


Orientations 

In order to produce a range of gravity effects while in 
terrestrial gravity, the evaporator was configured in five 
orientations: horizontal heater up (HU), vertical upflow (VU), 
horizontal heater down (HD), vertical downflow (VD), and 
sideways (SW), as shown graphically in Figure 8. 
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Figure 8: Evaporator orientations 


The appropriate length scales used in the calculation of 
the Edtvés number remain a subject of debate, as discussed in 
the "Microscale Definitions" section. Table 1 lists the Eétvés 
number for each orientation using three length scales: 

1. Conventional definition of the hydraulic diameter, 
squared; 

2. Cross sectional area of the microgap channel, H -W; 
and 

3. Channel width, W, for the surface tension term and 
length in the direction of the gravity vector (i.e., H for 
the HU and HD configurations, L for the VU and VD 
configurations, and W for the SW configuration) for the 
gravity term. 

The calculations revealed that: 

1. Eo is constant across all five orientations when calcula- 
tion methods | and 2 are employed; 

2. Calculation method 3 shows a strong dependence on 
orientation, producing results that vary by a factor of 60; 
and 


3. Calculation methods 2 and 3 produce identical results 
for the horizontal configurations, which results from the 
microgap height being in the direction of the gravity 
vector in those configurations. 

4. Based on the various proposed Eétvés number criteria 
for the transition to microscale, the microgap may or 
may not behave as a microchannel and therefore, may or 
may not provide gravity-insensitive flow boiling per- 
formance. 


Table 1: Eé6tvés number varies by factor of 900 for fixed 
geometry depending on calculation method 


Orientation E6tvés number 
HU Ap:g'H:-W 
Oo 
HD = 39 
VU Ap-g:Dp Ap-g:'H-W Ap-g:'L-W 
Oo oO Oo 
VD = 0.25 = 3.9 = 227 
Ap:g:'W:-W 
SW o 
= 233 


Characterization of Heat Losses to Ambient 

Prior to performing two-phase flow studies, the test facili- 
ty was characterized to assess the thermal losses to ambient. 
Two thermal resistance paths control the temperature of the 
TTC. Both paths begin at the resistors located on the frontside 
of the TTC. The upward path is a series of conduction through 
the silicon TTC and convection into the fluid flow, with 
thermal resistances defined as 


Rup = Reona,si + Reonv (14) 
tsi 
R , = ————_ 
cond,Si Kj z Achip (15) 
1 
| ee 
soit reas (16) 


where R,, is the thermal resistance of the upward path, 
Reona,si 18 the thermal resistance due to conduction through 
the silicon, ts; is the thickness of the silicon, Ks; is the thermal 
conductivity of the silicon, A;npip is the area of the TTC 
exposed to the fluid flow, Rony is the thermal resistance due 
to convection to the fluid, and HTC is the heat transfer 
coefficient. The downward path is a series of conduction 
through the solder and underfill, conduction through the 
printed circuit board, and conduction through the stainless 
steel flow enclosure. The complexity of the downward path 
makes the thermal resistance difficult to estimate analytically. 
In order to estimate the resistance of the downward path, the 
flow enclosure was evacuated to vacuum, heat was added to 
the TTC, and the temperature difference between the TTC and 
flow enclosure was recorded after steady state was achieved. 
The temperature difference was measured as the TTC power 
was increased to 5 W in | W increments and while the power 
was subsequently decreased in 1 W increments. 


A linear regression was applied to each set of temperature 
versus power data. The slope of each regression represents the 
thermal resistance between the diode and stainless steel 
enclosure. Table 2 provides the thermal resistance of the 
downward path for each diode. Diode 6, located at the center 
of the TTC, has the highest thermal resistance and diode 10, 
located immediately adjacent to the stainless steel housing, 
has the lowest thermal resistance. Diode 2 is representative of 
the diodes that fall between the maximum and minimum. 


Table 2: Downward thermal resistance for each diode 


Thermal Correlation 
Diode Resistance Coefficient 
(K/W) Squared 

1 13.42 0.9998 
2: 13.89 0.9998 
3 14.05 0.9998 
4 13.62 0.9998 
5 14.27 0.9997 
6 14.68 0.9997 
7 14.53 0.9998 
8 13.81 0.9998 
9 13.57 0.9998 
10 13.06 0.9998 


Figure 9 shows the fraction of the TTC heater power that 
enters the fluid and the fraction of the convective resistance to 
the total upward resistance as a function of heat transfer 
coefficient. The calculations assume a uniform heat transfer 
coefficient across the TTC and a constant silicon thermal 
conductivity of 113.7 W/m-K at 90 °C. The variation in 
downward thermal resistance among the diodes has a minimal 
effect on the fraction of the heat entering the fluid. For diode 
2, which represents an average diode on the TTC, 80% of the 
heat enters the fluid for heat transfer coefficients above 1,800 
W/m?-K, 90% of the heat enters the fluid for heat transfer 
coefficients above 4,100 W/m7-K, and 95% of the heat enters 
the fluid for heat transfer coefficients above 8,900 W/m’-K. 
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Figure 9: More than 80% of the heat enters the fluid for most of 
the operational space 


Based on the characterization study, the heat lost to am- 
bient (q,9s;) during two-phase testing was estimated using 


n [Gi—ts) 


Raowni 


dioss = » a ae 


i=1 


(17) 


where T; is the temperature of i-th diode on the TTC, T,, is the 
exterior temperature of the stainless steel fluid enclosure, 
Raown, 18 the downward thermal resistance of the i-th diode, 
and n is the number of diodes. This calculation method 
introduces a small error, as the diodes are not uniformly 
distributed across the area of the TTC. 
The temperature drop due to conduction through the TTC 
(AT ona) Was estimated for each diode using 
AT ona,i = Reona,si,i * (4 — oss) (18) 
where q is the total TTC power. This calculation assumes 
one-dimensional conduction. The single- and two-phase heat 
transfer coefficients (HTC) were calculated using 


q — oss 


HTC = 
Achip . (T; _ AT condi 


7) (19) 


where T; is the local fluid temperature. For single-phase flow, 
the local temperature of the fluid was assumed to increase 
linearly from the microgap inlet to outlet. For two-phase flow, 
the local temperature of the fluid was calculated based on the 
estimated local saturation pressure. It was assumed that the 
fluid pressure decreased linearly from inlet to outlet. For the 
highest flow rate and highest heat flux, the maximum two- 
phase pressure drop was 8.96 kPa, which corresponds to a 
drop in the local saturation temperature of 2.7 K. 


Operating Sequence 
The following sequence was used while collecting the 
test data. The data acquisition rate was 25 Hz. 
1. Enable reservoir heater and allow temperature and 
pressure to stabilize 
2. Enable recirculating chiller 
3. Apply power to pump to circulate working fluid at 
desired flow rate 
4. Apply power to preheater to achieve desired inlet fluid 
temperature 
5. Adjust the pump control voltage, preheater power, and 
reservoir set point to achieve desired conditions at the 
evaporator inlet 
6. Record data for 10 seconds when the flow loop has 
stabilized with the desired conditions at the evaporator 
inlet 
7. Apply 1 W/cm’ to the TTC, allow the TTC temperatures 
to stabilize, and record data for 10 seconds 
8. Increase the power to the TTC in 1 W/cm? increments, 
recording data for 10 seconds at each increment, until 
vapor is generated on the TTC 
9. After the onset of boiling, adjust the reservoir pressure, 
preheater power, and pump control voltage to maintain 
constant saturation pressure, degree of subcooling, and 
flow rate at the evaporator inlet 
10. After the loop has stabilized, record data for 10 seconds 
and capture photos of the two-phase flow 


11.Increase the heat flux by 1 W/cm”, adjust the reservoir 
temperature, preheater power, and pump control voltage 
to maintain constant saturation pressure, degree of 
subcooling, and flow rate at the evaporator inlet 

12. Repeat step 10 

13.Repeat steps 11 and 12 until the over temperature pro- 
tection limit is hit, increasing the heat flux increment up 
to 4 W/cm’ based on experience (however, the heat flux 
increment should be 1 W/cm’ as the maximum heat flux 
is approached) 

14.After the TTC temperature limit is reached, return the 
flow loop to the initial pump control voltage, preheater 
power, and reservoir set point, allow the flow loop to 
stabilize, and record data for 10 seconds 


Results and Analysis 

Single-phase and two-phase data points were collected at 
five orientations and three mass fluxes ranging from 349.1 to 
853.5 kg/m?-s. For each combination of mass flux and 
orientation, the heater power was increased until the 
overtemperature limit of 120 °C was reached. 

Figure 10 shows the critical heat flux (CHF) achieved at 
three mass flux increments with the evaporator in five 
orientations. CHF increased with increasing mass flux for all 
orientations, as expected. The horizontal heater up, vertical 
upflow, and sideways configurations produced CHF values 
that were 4.6 to 5.3% greater than those achieved in the 
horizontal heater down and vertical downflow configurations, 
which may be attributable to the differing role of gravity 
across orientations. In the horizontal heater up and vertical 
upflow configurations, gravity assists with vapor removal by 
accelerating it along or above the heated surface. In the 
horizontal heater down configuration, gravity pushes the 
vapor toward the heated surface and in the vertical downflow 
configuration, gravity pushes the vapor against the bulk fluid 
flow. It is also possible that the variation in the mass flux 
across orientations contributed to the variation in CHF, with 
the lower mass flux values in the horizontal heater down and 
vertical downflow configurations contributing to the lower 
values of CHF. Across the five orientations, the maximum to 
minimum variation in mass flux was 6.2% of the average at 
the lowest mass flux (361.8 kg/m’-s), 3.2% of the average at 
the intermediate mass flux (593.3 kg/m7-s), and 3.6% of the 
average at the highest mass flux (817.9 kg/m”-s). 

Figure 11 shows the single- and two-phase heat transfer 
coefficients at the lowest mass flux (left column) and highest 
mass flux (right column) for three diodes located along the 
centerline of the microgap. The top row shows the heat 
transfer coefficients (HTCs) for the diode located near the 
microgap inlet, the middle row shows the HTCs for the diode 
located in the middle of the microgap, and the bottom row 
shows the HTCs for the diode located near the microgap 
outlet. The distance between the first and second and second 
and third diodes is 6.35 mm (for a total length of 12.70 mm). 
As expected, the single-phase HTCs are much lower than the 
two-phase HTCs. At the onset of nucleate boiling, the heat 
transfer coefficient more than doubles for both mass fluxes. 
At the low mass flux, the vertical upflow orientation required 
the least superheat to initiate boiling and at the high mass flux, 
the horizontal down orientation required the least superheat to 


initiate boiling. The lack of clear trends in the required 
superheat to initiate boiling suggests that orientation is not the 
dominant mechanism affecting the onset of nucleate boiling in 
the present study. 
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Figure 10: Flow boiling critical heat flux varies negligibly with 
orientation across three mass fluxes 


Despite some variation in the onset of nucleate boiling, 
the heat transfer coefficients for each diode at each mass flux 
generally agree well across the five orientations studied, as 
evidenced by the excellent overlap of the five curves in each 
plot over much of the heat flux range. The biggest discrepan- 
cy appears to be for the diode near the channel inlet for the 
sideways configuration at the high mass flux, with slightly 
lower heat transfer coefficients than the other configurations. 
It is possible that this deviation resulted from a discrepancy in 
the inlet condition across orientations, as the heat transfer 
coefficients calculated based on the downstream diodes agree 
well with those for the other configurations and the CHF 
values for the sideways configuration were generally close to 
those for the preferable horizontal heater up and vertical 
upflow configurations. 

Regarding the effect of axial position on heat transfer 
coefficient, both mass fluxes showed the highest heat transfer 
coefficients at the microgap inlet across much of the range of 
heat fluxes. In some cases, at low heat fluxes, vapor bubbles 
did not nucleate until some distance downstream of the 
microgap inlet; in those cases, the downstream heat transfer 
coefficients are higher than those at the microgap inlet until 
the heat flux increases to the point that vapor is generated 
closer to the microgap inlet. For all mass fluxes, the highest 
overall heat transfer coefficients are achieved near the 
microgap inlet just before CHF. At each mass flux, as CHF is 
approached, the heat transfer coefficient near the microgap 
outlet begins to depart from the straight line expected during 
nucleate boiling. No such change was seen in the upstream 
heat transfer coefficients. The deterioration could be indica- 
tive of the transition to the intermittent flow regime, which is 
characterized by lower heat transfer coefficients than those 
achieved during bubbly flow. 
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Figure 11: Single- and two-phase heat transfer coefficients for HFE7100 in a 218 ym tall by 13.0 mm wide microgap as a function of 
evaporator orientation, mass flux, and sensor position 


Conclusions and Future Work 

A study has been performed to assess the role of gravity 
on two-phase flows in miniature and microscale rectangular 
ducts. A comprehensive literature review revealed that many 
of the criteria proposed for the transition to microscale two- 
phase flow, which is also the point at which the role of gravity 
is negligible, could be simplified to constant values of the 
Eétvés number. The wide variation among the proposed 
criteria may indicate that the transition cannot be captured by 
surface tension and gravity terms alone. Recent studies have 
incorporated the contact angle, mass flux, and fluid viscosity 
with some success. 

The orientation-dependence of the two-phase thermofluid 
behavior of HFE7100 was studied at three mass fluxes in a 
microgap cooler with a height of 219 um and width of 13.0 
mm. The results revealed that gravity played a negligible role 
in the CHF and heat transfer coefficients achieved along the 
axial length of the microgap cooler. These results further call 
into question the use of the Eétvés number for assessing the 
transition to microscale behavior since — depending on the 
formulation used — the range of Edtvés number for the present 
system ranged from 0.25 to 233 across the five orientations 
studied. Additional work is required to assess the utility of the 
E6tvés number for predicting the microscale transition and/or 
to determine the appropriate length scales for calculating the 
E6tvés number. 

Future work, using the existing test facility, will include 
microgap heights from 100 to 1000 um and mass fluxes from 
100 to 2000 kg/m’-s to better assess the effects of microgap 
height and fluid velocity on the orientation-dependence of the 
two-phase thermofluid behavior. Additional studies will be 
performed to assess other parameters, such as channel aspect 
ratio, heat flux, vapor quality, and/or flow regime. The role of 
flow regimes is of particular interest, as such a dependence 
could imply that each regime has a unique criterion and 
therefore knowledge of the flow regime would improve the 
predictive accuracy of the criteria. Bar-Cohen and Rahim 
previously found that a regime-informed approach improved 
the predictive accuracy of correlations for two-phase heat 
transfer coefficients [35]. It is also possible that gravity is 
important to only certain regimes and their transitions. For 
example, Taitel [36] found that in the case of the transitions to 
annular flow and dispersed bubble flow, the effect of gravity 
was negligible. 

Ultimately, tests in a microgravity environment will be 
required to validate the ground-based test results. Microgravi- 
ty validation would then allow the criteria that result from the 
research program to inform the design of gravity-insensitive 
flow boiling systems for future space missions. 


Nomenclature 
A Area (m”) 
Bo _ Bond number (-) 
Co Confinement number (-) 
D Tube diameter (m) 
Dn Hydraulic diameter (m) 
Eo _ E6tvés number (-) 
F Force (N) 
G Mass flux (kg/m7-s) 
H Channel height (m) 


HD _ Horizontal heater down 
HU _ Horizontal heater up 
Heat transfer coefficient (W/m’-K) 
L Channel length (m) 
La __ Laplace constant (-) 

P Channel wetted perimeter (m) 

R Thermal resistance (K/W) 
Re __ Reynolds number (-) 
Sideways 

T Temperature (°C) 
VD _ Vertical downflow 
VU _ Vertical upflow 

Ww Channel width (m) 

g Acceleration due to gravity (m/s”) 
k Thermal conductivity (W/m-K) 
q 

r 

t 


Heat (W) 
Tube radius (m) 
Thickness (m) 
Subscripts 
Si Silicon 
c Cross-sectional 
cond Conduction 
f Fluid 
g Gas 
l Liquid 
ss Stainless steel 
v Vapor 
Greek Letters 
a Ratio of gravitational forces to surface tension 
forces (-) 
6 Contact angle (°) 
L Dynamic viscosity (kg/m-s) 
p Density (kg/m?) 
0 Surface tension (N/m) 
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